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Abstract 

I define a model of three-dimensional simplicial gravity using an extended 
^ ' ensemble of triangulations where, in addition to the usual combinatorial tri- 

, angulations, I allow degenerate triangulations, i.e. triangulations with distinct 

' simplexes defined by the same set of vertexes. I demonstrate, using numeri- 

cal simulations, that allowing this type of degeneracy substantially reduces the 
geometric finite-size effects, especially in the crumpled phase of the model, in 
QQ , other respect the phase structure of the model is not affected. 
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1 Introduction 

Among the different attempts to construct a viable theory of quantum gravity, 
Euclidean quantum gravity is one of few that allows for non-perturbative studies 
r> I using numerical methods. In the Euclidean approach one works with a functional 

' integral over Euclidean metrics g weighted by the Einstein-Hilbert action: 

SE{[9ab]) = j d^^y^{R-2A), (1) 

where G is the Newton's gravitational constant, A the cosmological constant and R 
the scalar curvature of the metric gab on the space-time manifold M. 

Simplicial gravity, or dynamical triangulations, is a regularization of this 
model where the integration over metrics is replaces by all possible gluings of D- 
simplexes into closed simplicial manifolds [Q, |2[ . In three dimensions the regularized 
partition function becomes 

Zi^.,^) = Y: ^e-^^3 + KiVo (2) 
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where T denotes an appropriate ensemble of three-dimensional simplicial manifolds, 
Ct is the symmetry factor of a triangulation T, iVj is the number of z-simplexes in 
the triangulation, and and k are coupling constants related to the cosmological 
and the inverse Newton's constant, respectively. 

As it stands the model Eq. (|2|) needs not to be well defined for any value 
of the couplings; unless the sum over manifolds is restricted in some way it is not 
convergent |^]. The same is true in two dimensions where, without the restriction to 
triangulations of fixed topology, the number of triangulations is not exponentially 
bounded as a function of volume — a necessary condition for convergence of the 
sum. In three dimensions the situation is even worse, there is no proof that even 
restricted to fixed topology the sum is convergent. The existence of an exponential 
bound is, however, strongly supported by numerical simulations for an ensemble of 
spherical combinatorial triangulations |3|, Q . 

The model Eq. (|2|) has been studied extensively using numerical simulations 
0' 3' IB ^ discontinuous phase transition at a value of the inverse New- 

ton's constant k ~ 3.8, separating a strong-coupling (small k) crumpled phase from 
a weak-coupling branched polymer phase. The crumpled phase is dominated by 
triangulations of almost zero extension; characterized by few vertexes connected to 
a large portion of the triangulation and has a internal fractal dimension dn = oo. 
This excludes any sensible continuum limit in this phase as the distance between 
two points will always stay at the Planck scale. The branched polymer phase is, 
on the other hand, dominated by essentially "one-dimensional" triangulations — 
bubbles glued together along small necks into a tree-like structure — with fractal 
dimension two. 

Little attention has been paid so far to the ensemble of triangulations included 
in the partition function Eq. (^), although in two dimensions it is well known that 
while different ensembles yield the same continuum theory 0] the finite-size effects 
depend strongly the type of triangulations used. It has been shown that the less 
restricted the triangulations are — the larger the ensemble of triangulations is — 
the smaller the finite-size effects are ^. This result can be understood intuitively 
as, for a given volume, with a larger triangulation-space it is easier to approximate 
a particular fractal structure. This is especially important as in simulating models 
of dynamical triangulations the geometric finite-size effects usually dominate. 

Simulations of simplicial gravity in dimensions larger than two have until now 
used exclusively combinatorial triangulations. In a combinatorial triangulation each 
D-simplex is uniquely defined by a set of {D + 1) vertexes, i.e. it is combinatorially 
unique. In this paper I explore a different ensemble of triangulations, degenerate 
triangulations, defined by relaxing this restriction and allow distinct simplexes de- 
fined by the same set of vertexes. This allows, for example, two vertexes connected 
by more than on edge and simplexes that have more than one face in common. 
I retain, on the other hand, the restriction that for every simplex all its {D + 1) 
vertexes should be different, i.e. I exclude degenerate simplexes. As the ensemble 
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of degenerate triangulation includes combinatorial triangulations as a subclass, it is 
obviously larger and as I will demonstrate in this paper this leads to smaller finite- 
size effects, just as in two dimensions. Apart from this, using degenerate rather 
than combinatorial triangulations does not seem to change the phase structure of 
the model Eq. (^); as the coupling constant k is varied I observer a discontinuous 
phase transition separating a crumpled phase from a branched polymer phase. 

The paper is organized a follows: In Section 2 I define the different ensembles 
of three-dimensional triangulations to be considered and discuss some practical as- 
pects such as ergodicity of the algorithms used in the simulations and the possible 
problem with pseudo- manifolds. In Section 3 I describe the numerical setup and 
compare the efficiency of the simulations and the finite-size effects for the ensem- 
bles of degenerate versus combinatorial triangulations. In Section 4 I investigate 
the phase structure of the model, identify a phase transition and show that it is 
discontinuous, and explore the different phases. Finally, in Section 5 I discuss the 
practical consequences of this for further simulations of simplicial gravity. 

2 Degenerate triangulations 

2.1 Different ensembles of triangulations 

The triangulations included in the partition function Eq. (||) are constructed by 
successively gluing equilateral 3-simplexes (tetrahedra) together along their two- 
dimensional faces (triangles) into a closed three-dimensional simplicial manifold |l|, 
^. Pseudo- manifolds are eliminated from the ensemble by the restriction that the 
neighborhood of every vertex in the triangulation should be homeomorphic to a 
sphere. To each tetrahedra there is associated a set of four vertexes, or 0-simplexes; 
the order of a vertex (or more generally of a sub-simplex) is defined as the number 
of tetrahedra containing that vertex. A dual graph to a triangulation is defined by 
placing a vertex in the center of each tetrahedra and connecting vertexes in adjacent 
tetrahedra together. For a three-dimensional triangulation the dual is a (/)^-graph. 
A triangulation is said to be degenerate if it contains two distinct 3-simplexes, or 
sub-simplexes, which are combinatorially identical, i.e. defined by the same set of 
vertexes. A simplex is said to be degenerate if it is defined by a set of vertexes 
including the same vertex more than once. 

By imposing additional restrictions on how the simplexes may be glued to- 
gether, it is possible to define different ensembles of triangulations. In this paper 
I will consider two such ensembles; combinatorial and degenerate triangulations. 
Before discussing these different ensembles in three dimensions, let's start by de- 
scribing the corresponding triangulations in the simpler case of two dimensions. 
Given the definitions above we distinguish among three types of two-dimensional 
triangulations: 
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Figure 1: An example of degenerate two-dimensional triangulations. The dashed lines 
indicate the corresponding dual graphs. In (a) the triangles ta and th have two edges in 
common, {vi,V2} and {v2,V3} (a self-energy diagram); this type of degeneracy is allowed 
in restricted degenerate triangulations. In (6) the triangle ta = {vi,vi,V2} is degenerate 
(corresponds to a tadpole); this is allowed in maximally degenerate triangulations. 



(a) Combinatorial triangulations : Tc 

In a combinatorial triangulation every triangle, and every edge, is combinatori- 
ally distinct, i.e. it is uniquely defined by a set of distinct vertexes. From this it 
follows that no two triangles have more than one edge in common and a triangle 
cannot be its own neighbor. The dual graph to a combinatorial triangulation is 
a i;^^-graph, excluding both tadpoles and self-energy diagrams. 

(b) Restricted degenerate triangulations : Td^ 

By allowing distinct triangles (or edges) that are combinatorially equivalent, i.e. 
defined by the same set of three (two) vertexes, while still excluding degenerate 
triangles, we define a restricted ensemble of degenerate triangulations, Td^. 
This type of triangulations can have vertexes connected by more than one edge, 
and triangles that have two edges in common. In this case self-energy diagrams 
are allowed in the dual (/>^-graph. 

(c) Maximally degenerate triangulations : T^j^^ 

Finally, by also allowing degenerate triangles we define an ensemble of maxi- 
mally degenerate triangulations, Tq^. This ensemble allows self- loops — edges 
starting and ending at the same vertex — equivalently, both tadpoles and self- 
energy diagrams are included in the dual </>^-graph. 

Clearly Tq C Td^ C Tb^,^ . I show examples of those different types of degeneracy in 
Figure As models of pure two-dimensional simplicial gravity defined with those 
different ensembles can been solved analytically as matrix models, it is known that 
in all cases they define the same continuum theory Q. This is, of course, what 
one would expect based on universality; the different ensembles simply amount to 
different discretization of the manifolds and the details of the discretization should be 
irrelevant in the continuum limit. This is true even if the triangulations are restricted 
even further and only minimal curvature fluctuations of ±1 are allowed (including 
only vertexes of order 5, 6 and 7) Q. But, as mentioned in the Introduction, the 
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Figure 2: Two examples of degeneracy allowed in three-dimensional degenerate triangula- 
tions. Figures (02) and (62) show the corresponding dual (/(^-graphs. In both examples there 
are two distinct but combinatorially equivalent tetrahedra Sa and s;,, defined by the same 
set of vertexes {wi, V2-, t'a, ^"4}- In (ai) the tetrahedra have two triangles in common, whereas 
in (61 ) they share three triangles, in which case the dual graph is a self-energy diagram. 

finite-size effects do depend on the ensemble used, in particular they increase as 
more restrictions are placed on the triangulations [^, ^. 

Analogous to two dimensions, by imposing the appropriate restrictions in three 
dimensions I defined the corresponding three different ensembles of triangulations. 
By allowing distinct but combinatorially equivalent non-degenerate simplexes, I de- 
fine an ensemble of restricted degenerate three-dimensional triangulations. Com- 
pared to combinatorial triangulations, in this ensemble two tetrahedra may be glued 
together along more than one of their faces. I show examples of this in Figure ^. 

Relaxing these restrictions further and allow degeneracy within the simplexes 
themselves defines an ensemble of maximally degenerate triangulations|^, T^j^j This 
allows, for example, self-loops both in the triangulation and in its dual graph (a 
simplex that is its own neighbor). Although this type of degenerate triangulations 
is a priori well defined, my numerical investigation indicates that this ensemble does 
not lead to a well defined statistical model. I tried to simulate the model Eq. (|2|) 
using maximally degenerate triangulations; however, in quasi-canonical simulations 
(see Section 3) it was not possible to tune the cosmological constant /i in such a 
way that the volume fluctuated around a fixed target volume N^. Regardless of the 
value of n used, in the course of the simulations the volume eventually exploded. 

^Note that in contrast to two (or more generally even) dimensions it is not possible to construct a 
three-dimensional triangulation with tadpoles in the corresponding dual <j!>^-graph. For a part of the 
dual graph to be connected to the rest by a single link, the corresponding part of the triangulation 
would have to be enclosed by a (closed) boundary consisting of a single triangle. This is, however, 
impossible. 
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Figure 3: An example of a degenerate simplex, Sa = {vi,vi,V2,V3}, allowed in a three 
dimensional maximally degenerate triangulation. The edge Ca — is a self- loop. 




Figure 4: The (p, 9)-moves in three dimensions: (a) Inserting a vertex and its inverse, 
deleting a vertex. (6) Replacing a triangle by an edge and vice verse. 



This behavior suggests that the numbers of maximally degenerate three-dimensional 
triangulations grows faster than exponentially with the volume, in which case the 
partition function Eq. (Q) does not converge. In the rest of this paper I only consider 
the ensemble Tdri omitting the specification restricted when referring to it. 



2.2 Ergodicity in the updating procedure 

Having defined degenerate triangulations, a necessary prerequisite for using them in 
numerical simulations of simplicial gravity is that there exist a set of ergodic geo- 
metric moves which allow us to explore the triangulation-space. For combinatorial 
triangulations this is provide by the {p, g)-moves, a variant of the Alexander moves 



|10|. In three dimensions the (p, g)-moves consist of either inserting or removing 
a vertex of order four, or replacing a triangle by an edge of order three and vice 
verse. These moves are depicted in Figure To apply the same set of moves in 
simulations with degenerate triangulations I have to show that they are also er- 
godic when applied to the ensemble Tp. To do so it is sufficient to show that every 
degenerate triangulation can be changed into a combinatorial triangulation using 
a finite sequence of the {p, g)-moves. I distinguish among three different types of 
degeneracy in restricted degenerate triangulations: combinatorially equivalent tetra- 
hedra, triangles and edges, all of which can be eliminate following the appropriate 
procedures: 
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(a) Equivalent tetrahedra: 

In the case of two, or more, tetrahedra defined by the same set of vertexes it is 
sufficient to insert a vertex (applying move (1,4)) into each of those tetrahedra, 
this will create a set of combinatorially distinct 3-simplexes. 

(b) Equivalent triangles: 

Given two triangles defined by identical set of vertexes I can apply move (2, 3) 
to one of them, i.e. replace that triangle with a new edge. To ensure that the 
new edge does not result in multiply connected vertexes, and to make sure the 
move is not rejected by a geometric constraint, prior to using move (2,3) I can 
apply move (1,4) to both of the tetrahedra containing the triangle in question 
— the new edge now joins two previously unconnected vertexes. 

(c) Equivalent edges: 

It is slightly more complicated to eliminate a par of distinct edges connecting 
the same two vertexes. To remove an edge I must apply move (2, 3), i.e. replace 
it with a triangle. This is, however, only possible if the corresponding edge is 
of order three. To reduce the order of the edge I apply move (2,3) repeatedly 
to the neighboring triangles. It is possible that some of those triangles could be 
combinatorially identical; this I eliminate by first applying procedure (b). The 
remaining triangles can then safely be replaced until the order of the edge is 
three. 

An additional complication that could arise when restrictions on the triangu- 
lations is relaxed is the appearance of pseudo-manifolds In general when gluing 
i)-simplexes together into a closed manifold one does not get a simplicial manifold 
but a pseudo- manifold. A simplicial manifold is defined by the additional require- 
ment that the neighborhood of each vertex is homeomorphic to the D-dimensional 
ball Bd. As this requirement is trivially satisfied for Z) = 2, all degenerate trian- 
gulations in two dimensions are automatically simplicial manifolds. It is a priori 
not clear that in simulating degenerate triangulations in three dimensions, pseudo- 
manifolds could not be generated — something we would like to avoid. 

I observe, however, that in order to generate a pseudo-manifold the {p, q)- 
moves must change the topology of the dual complex to a given vertex vq in the 
triangulation. By the dual complex I mean the two-dimensional triangulation con- 
structed by connecting by edges the centers of all adjacent tetrahedra containing 
vq. But as the moves only change the triangulation, hence the dual complex, locally 
they will not change its topology. More precisely, when applying a {p, q)-move to 
a three-dimensional triangulation, degenerate or not, the change in the local neigh- 
borhood of each of the vertexes involved in the move can be given in terms of the 
corresponding two-dimensional (p, g)-moves (i.e. inserting/deleting a vertex or flip- 
ping an edge) applied to its dual complex. This ensures that the topology of the 
dual complex is not changed when the move is applied and pseudo-manifolds are 
not generated. 
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3 Simulating degenerate triangulations 



To study the model Eq. (^) I now turn to numerical simulations. As the {p, g)"moves 
in three dimensions do not preserve the volume of the triangulation the simulations 
cannot be restricted to the canonical ensemble, contrary to what is possible in 
two dimensions. In addition, fluctuations in the volume are necessary to maintain 
ergodicity in the updating procedure. In practice though I perform quasi-canonical 
Monte Carlo simulations, with almost fixed N3, simulating the partition function 



Z{^,.;N,) = ^e-/^^3-7(iV3-iV3)^f^^^(^), (3) 



tgT{N3) 



Here Q]\i^{k.) is the canonical partition function and the quadratic potential term 
added to the action ensures, for an appropriate choice of the parameter 7, the 
necessary fluctuations around a target volume N^. In simulating degenerate trian- 
gulations I found 7 ~ 0.001 to be an appropriate choice. In all cases I simulate an 
ensemble of spherical triangulations. 

From a practical point of view simulating degenerate triangulations is actually 
simpler than simulating combinatorial triangulations. For the latter the most time- 
consuming part of the update is the manifold check, i.e. to verify that there exist no 
sub-simplex in the triangulation combinatorially equivalent to the one created by the 
move [11]. In three dimensions this check is done both when a triangle is replaced 
by an edge, and vice verse. This check is not needed for degenerate triangulation; 
the only geometric check that remains is in replacing a triangle by an edge the two 
tetrahedra involved must only have one triangle in common. As this check is local 
the computational effort needed to update a degenerate triangulation is less than 
for a combinatorial one. 

This does not automatically imply that simulations with degenerate triangu- 
lations are more efflcient. As the ensemble of degenerate triangulations is much 
larger the corresponding critical cosmological constant /Xc is also larger than for the 
ensemble of combinatorial triangulations. And as fic enters the detail balance con- 
dition used in the Metropolis test this in effect reduces the acceptance rate in the 
updating procedure. To compare the efficiency of simulating the different ensembles 
I must estimate how fast, measured in CPU-time, independent configurations are 
generated. I have done this for both ensembles in the crumpled phase, k = 0. I 
caution, however, that the comparative efficiency will depend on several factors: the 
particular implementation of the algorithm, how efficiently the manifold checks are 
executed, the probability of choosing the different moves, and where in the phase- 
space the simulations are carried out. In the comparison I use the same program 
for both types of triangulations, turning on and off the manifold checks, and choose 
the different moves with equal probability. I measure the CPU-times tg (in ms) it 
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Figure 5: The effective critical cosmological constant ndN^), both for an ensemble of 
degenerate (circles) and combinatorial (diamonds) triangulations. The solid lines are fit to 
the form: ndN^) — p,c + P/N^ ■ Insert: The relative deviation from the infinite volume 
value: Afic = \p.c - ^c(A^3)|/Mc. 

takes to perform one "sweep" through the triangulations; a sweep is defined as 
attempted local moves, and the auto-correlation times r in units of sweeps, mea- 
sured for the time-series of the energy-density ng = Nq/N^. Combining tg and r, I 
get an estimate of how fast independent configurations are generated: 



For small triangulations, ~ 800, the effective auto-correlations are comparable for 
the two ensembles, whereas for larger volumes simulating degenerate triangulations 
is substantially faster. 

The real benefit of using degenerate triangulations comes from the reduction 
in finite-size effects. This reduction is especially pronounced in the crumpled phase 
where the internal geometry is dominated by collapsed manifolds. To demonstrate 
this I measured the volume dependence of the effective critical cosmological constant 
fij?{Ns) for K = and for triangulations of volume < 25000. This I show in 
Figure |. For comparison I include in the figure the corresponding values jjL^{N^) 
for combinatorial triangulations. For both ensembles f^dN^) approaches a finite 
value flc as the volume diverges; this in turn implies that the canonical partition 
function fi^Vg is exponentially bounded as a function of volume — an necessary 
condition for the partition function Eq. (^) to be convergent. I have fitted the data 
to an assumption of a power-law convergence to an exponential bound (including 




0.43 N^^-^'^ combinatorial, 
230 N^^-^"^ degenerate. 



(5) 
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volumes > 400); this yields: 

2.073(21) - 3.67(44) Nf^-^^^'^\ = 4.0 combinatorial, 
/ic(A^3) = { (6) 
2.447(3) - 2.97(29) iVg^^-^^^^S = 0.8 degenerate. 

For both ensembles the quality of the fit (the x^/d.o.f) is acceptable. This can 
be compared to a fit to a super-exponential growth of 0,]\f^, compatible with the 
assumption /ic(-^3) ~ a + piogN^, which is definitely ruled out by a very large 
value. 

There is, however, a marked difference in how fast fic approaches its infinite 
volume value for those two ensembles. This is shown in the insert in Figure ^ 
where I plot Afic = \p-c — /^c(-^3)|//^c. Taken at face value, this suggests a reduc- 
tion of finite-size effects by two orders of magnitude when using degenerate rather 
than combinatorial triangulations. This huge reduction is though presumable only 
achieved in the crumpled phase where collapsed manifolds dominate the partition 
function and the geometric finite-size effects are most pronounced. 



4 The phase structure 

4.1 Existence of a phase transition 

I now turn to the existence of a phase transition in the model as the inverse New- 
ton's constant k is varied. For combinatorial triangulations it is well established by 
numerical simulations that the model Eq. (^) has a discontinuous phase transition 
at K ~ 3.8, separating a strong-coupling (small k) crumpled phase from a weak- 
coupling branched polymer phase j^, H] . As I will demonstrate in this section this 
is also the case for degenerate triangulations; moreover, the discontinuous nature of 
the transition is easily observed on triangulations of relatively modest size. 

I have simulated the model Eq. (^ for target volumes N3 = 200 to 6400, and 
for each volume I search for a phase transition in the fluctuations of the energy- 
density riQ = Nq/N^. a signal for a transition would be a peak in the specific 
heat: 

Cn, = ^ {{N^) - {No)') . (7) 

This is shown in Figure |6| for different target volumes. I did simulations at few 
values of k around the observed peaks, collecting between two and ten thousand 
independent measurements, and then used standard multi-histogram methods to 
interpolate between the measurements. The interpolations is shown as solid lines in 
Figure |6|, the dashed lines indicate the errors. 

To infer about the nature of the transition I look at the energy-density dis- 
tribution p(no). For sufficiently large target volume this distribution has a well 
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Figure 6: The measured specific heat Cn^ — the fluctuations in tlie number of vertexes iVg 
— as the inverse Newton's constant k, is varied. This is for target volume A^3 — 400 to 6400. 
The solid lines are interpolations between measurements, using multi-histogram methods, 
and the dashed lines indicate the error. 



resolved double-peak structure, indicating a discontinuous phase transition. This 
double-peak structure is observed already on volumes > 1600. I show an exam- 
ple of the measured energy distributions in Figure (thin lines) for the three largest 
volumes I studied, together with fits to a form composed of two Gaussian peaks 
(thick lines): 

p{no) = ai e-'^i (^1-^0)% ^2 e-^2 (62 - no)^ _ (g) 

Each distribution is normalized by the height of the peaks and the value of k is 
chosen so that the two peaks are of equal height. 

Additional evidence for a phase transition comes from other geometric observ- 
ables such as the maximal order of a vertex, qq. As discussed in the Introduction, 
the crumpled phase is dominated by triangulations containing few vertexes of very 
large order, whereas in a branched polymer phase the vertex orders are more or 
less equally distributed. This is reflected in a sudden drop in the maximal vertex 
order across the transition; this change leads to divergent peaks in both the fluctu- 
ations of qo, i.e. the susceptibility Xqo = -^3 ((^o) ~ (9o)^)i and its energy derivative 
Dqp = N3 {{qono) — ((i'o)(^o))- In fact, plots of both Xqo and Dg^ look more or less 
identical to Figure |6|. For all the observables the peak heights diverge, an example 
of this is shown in Figure ^ for the specific heat. I have estimated the corresponding 
scaling behavior, fitting the peak values to the form: Opeafe ~ A*""; the fitted scaling 
exponents are shown in Table 1. Compared to Xqo and Dgp, the fit of the specific 
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Figure 7: The energy- density distribution p{no) for N3 = 1600, 3200 and 64000. The thin 
Unas are the measured distributions, the thick hues are fits to two Gaussian peaks Eq. (H). 
Each distributions is normalized with the height of the peaks and the values of k chosen so 
the two peaks are of equal height. 



heat is rendered considerable more difficult by an additional constant contribution 
to the scaling behavior. Nevertheless, for the specific heat the peak height scales 
almost linearly for sufficiently large target volume as is expected for a discontinuous 
phase transition. 

To locate the infinite volume critical coupling Re I fit the pseudo-critical cou- 
plings KciN^), defined by the location of peaks in the different observables, to the 
expected finite-volume scaling behavior: 

, ,^ , const. , , 

Kc(A^3) = + (9) 

The fitted values of Kc and the exponent x are shown in Table 1. To asses the finite- 
size effects I repeated the fits with a different lower cut-off on the included target 
volume. This did not affect the estimate of Kc appreciatively; in all cases Kc ~ 3.8. 
Note however that, in contrast to what has recently been reported in simulations 
with combinatorial triangulations I observe no indication of the peak locations 
running to Kc = 00. For all the observables the pseudo-critical couplings converge 
nicely. This is shown in Figure |9[ A fit to an assumption of non-convergence, 
i^d^s) = a + blog N^, is ruled out by a x^/(d.o.f.) > 1000. For the scaling exponent 
X, I get values x ~ 1/3. For a discontinuous phase transition one would naively 
expect 6kc = |kc(-^3) — ~ L^'^^ , where L is some appropriate length scale in the 
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Figure 8: The volume scaling of tiie peak in the specific heat. The solid line is a fit to the 
scaling form: Cn^ ~ a + bN^. The value of the exponent a is shown in Table 1. 

system and dn is the (internal) fractal dimension (see e.g. Ref. fl^). This would 
imply X = 1 in Eq. (^). However, for a system with a fluctuating internal geometry 
like simplicial gravity, it is not clear what dimensionality should be associated to the 
system at a discontinuous phase transition where in one phase the internal fractal 
dimension is infinite, but two in the other. Hence it is not obvious how to interpret 
the observed value of the exponent x. 

4.2 The crumpled phase 

As already stated the strong-coupling phase of the model Eq. (|2|) is a crumpled 
phase, characterized by few vertexes of very large order. Qualitatively 1 observe 
this behavior regardless of the type of triangulations 1 simulated. Note, however, 
that there is a priori no reason to expect exactly the same behavior for the two 
ensembles in this phase; because of the collapsed nature of manifolds it is unlikely 
that any sensible continuum limit can be taken for k < Kc- Hence the details of the 
triangulations — the discretization — may be relevant in the infinite volume limit. 

There does indeed appear to be some slight difference in the internal geometry 
of the ensemble of degenerate triangulations, compared to the combinatorial ensem- 
ble, for K = 0. An example of this is the scaling behavior of the effective critical 
cosmological constant determined in Section 3. For combinatorial triangulations the 
exponent of the power-law correction is 7 ~ 0.3, whereas 7 ~ 0.5 for degenerate tri- 
angulations. Similar difference is observed in the volume scaling of the order of the 
most singular vertex. For three-dimensional combinatorial triangulations go scales 
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Table 1: (A) The volume scaling of the peaks in C 



Ns^ Xqo 



bN: 



3 , together with the quality of the fits. 



and Dqg , fitted to the form: 
For Xqo ^'^d T)qQ I use a — 0. Values are 
on the target volume include in the fit. (B) The 
corresponding values of Re, determined from fitting the pseudo-critical coupling constants 
Kci^N^) to the form Eq. (|^), together with the scaling exponent x. 
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Figure 9: Scaling of the pseudo-critical inverse Newton's constant Kc{Nz), defined by the 
location of peaks in Cn^, , Xqo ^^"^ ^qo respectively. The values are plotted versus the volume 
scaled with an exponent x — 1/3. 
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Table 2: The fractal dimension d^j for an ensemble of degenerate three-dimensional tri- 
angulations, determined by collapsing the simplex-simplex correlation function gjva {^) onto 
the scaling form Eq. (10). The values shown are for k = in the crumpled phase, and 
K = 5 and 8 in the branched polymer phase. The analysis is done including different target 
volumes N^^. Also included is the x^/(d-0.f.) value of the collapse. 





K = 
dH 


K = 5 
dH X^ 


K = 8 
dn X^ 


{200 — 400} 
{400 — 800} 
{800 — 1600} 
{1600 — 3200} 
{3200 — 6400} 

{400 — 6400} 


4.05(15) 55 
4.66(18) 62 
5.27(14) 50 
5.68(16) 32 
6.25(20) 30 

5.20(30) 84 


1.85(6) 0.64 
1.92(8) 0.84 
1.93(8) 0.77 
2.00(13) 0.45 
1.99(15) 0.45 

1.94(10) 1.1 


1.89(5) 0.55 
1.95(7) 0.55 
1.93(13) 0.32 
1.98(15) 0.58 
1.96(12) 0.50 

1.95(11) 0.6 



42f2) 

sub-linearly with volume: qo ~ N^' , for degenerate triangulations, on the other 
hand, the scaling exponent is larger; ~ A'^^'^'^'-^^ 

The collapsed nature of this phase is most easily demonstrated by measuring 
the fractal dimensions dn- To estimate dn I measure the two-point correlation 
function gN^{r), defined as the number of simplexes at a geodesic distance r from 
a marked simplex. The geodesic distance between two simplexes is defined as the 
shortest path traversed through the center of simplexes. Assuming that the only 
relevant length scale in the system is given by N^^^" , general scaling arguments |jl^ 
imply that 

gN,{r) ~ Nl-'/'" Fix), x = (10) 

^3 

The fractal dimension is measured by "collapsing" distributions g^^ir), correspond- 
ing to different volumes, onto a single scaling curve. I have done this for volume 
A^3 < 6400 and for k = 0; the value of dn, determined by minimizing the x^^value 
of the collapse, is shown in Table 2. The errors are estimated using jack- knifing. 
To demonstrate how the estimate of dff depends on the volume, I have done the 
analysis either including all the volumes, or using pairs of consecutive volumes A''3 
and 2A'^3. For n = the estimate of dn drifts systematically toward larger values 
as the volume is increased, in addition the fit to the scaling behavior Eq. (|To|) is 
very poor; the x^/(d.o.f.) >> 1. This is not unexpected, in a crumpled phase the 
internal fractal dimension is dn = oo, in which case the simplex-simplex correlation 
function is not expected to scale. 
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Figure 10: The scaled simplex-simplex correlation function gj^^{x) = gN3ix)/N^ 
versus the scaling variable x = r/N^^'^" . This is for degenerate three-dimensional triangu- 
lations of size N-^ = 400 to 6400 and k = 5. The measured value of dn = 1-94 has been used 
in the scaling. 



4.3 The branched polymer phase 

Finally I have also investigated the nature of the weak-coupling phase for degenerate 
triangulations at two values of the coupling constant: /t = 5 and 8. To establish 
that this is a branched polymer phase I measured both the fractal dimension dn 
and the string susceptibility exponent 7^. The estimate of dn is included Table 2. 
In contrast to the crumpled phase in this phase I get a consistent value dn ~ 2, 
as expected for branched polymers; for both values of k the quality of the scaling 
is very good, x^/(d-o.f.) < 1. An example of the scaling is shown in Figure |l^ for 
K = 5. 

The string susceptibility exponent 7^ defines the sub-leading correction to the 
large-volume behavior of the canonical partition function: ~ exp{ficN3}N^"~^ ■ 
An efficient method for determining 7^ in quasi-canonical simulations is provided by 
the size distribution of baby universes. A baby universe is a part of a triangulation 



connected to the rest via a small "neck", usually chosen as a minimal neck |15|. 
For a three-dimensional combinatorial triangulation a minimal neck consists of four 
triangles glued together along their edges into a closed surface (a tetrahedra). By 
cutting the triangulation up along this boundary it is divided into two parts; a 
"mother" of size A^3 — B and a "baby" of size B. The size distribution of baby 
universes, bc4^{B), is defined as all such partitions averaged over the ensemble. From 
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the asymptotic behavior of the canonical partition function it follows that 



The distribution bc^^iB) can be measured with high accuracy in quasi-canonical 
simulations and the exponent 7^ is extracted by a fit to Eq. (11). 



The above definition of a baby universe also applies to degenerate triangu- 
lations although in this case a minimal neck consist of only two triangles glued 
together along their three edges (see e.g. Figure 26). This defines a distribution of 
baby universes, b£)^{B)^ separated from the rest by two links in the dual graph. It is 
also possible to define a distribution of baby universes, boiiB), separated from the 
rest by four links in the dual graph, analogous to the definition for combinatorial 
triangulations. I have verified that the two distributions, boiiB) and boiiB), lead 
to the same estimate of 7^. 

I measured the distributions for k = 5 and on volumes A^3 = 200, 400 and 800, 
and estimated 7^, by a fit to Eq. (pT|). This gave the following estimates: 7^ = 0.47(3), 
0.48(3) and 0.44(4), for = 200, 400 and 800, respectively. This agrees reasonable 
with the expected value for branched polymers, 7^ = 1/2, especially given the small 
size of the triangulations used in the estimate. 



5 Discussion 

I have investigated a model of three-dimensional simplicial gravity defined with 
an ensemble of degenerate triangulations allowing distinct simplexes sharing the 
same set of vertexes. The motivation for using degenerate triangulations, instead 
of the combinatorial ones traditionally used in simplicial gravity for D > 2, is 
the experience from two-dimensional simplicial gravity which shows that using less 
restricted triangulations, i.e. a larger triangulations-space, substantially reduces the 
finite-size effects in numerical simulations. 

As in two dimensions using a different ensemble of triangulations, provided it 
yields a well defined statistical model, does not seem to change the phase structure 
of the model; I still observe a crumpled and a branched polymer phase, separated 
by a discontinuous phase transition. There is, however, evidence of less finite-size 
effects in the model, especially in the crumpled phase where the critical cosmologi- 
cal constant approaches much faster its infinite volume value. As a consequence, an 
exponential bound on the number of degenerate triangulations of a given volume is 
more easily established than for the corresponding ensemble of combinatorial trian- 
gulations. This reduction in finite-size effects in the crumpled phase is intuitively 
clear; with less restrictions on the connectivity of the simplexes the degenerate tri- 
angulations crumple more easily. In addition, the discontinuous nature of the phase 
transition is easily observed on triangulations of relatively modest size. 
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A natural extension of the work presented in this paper is to add a measure 
term [^], or alternatively matter fields [17|, to the action Eq. and to investigate 
if the phase structure of the model is modified in the same way as has recently 
been observed in simulations with combinatorial triangulations in both four jl^ 
and three dimensions . There a new phase — a crinkled phase — appears for a 
sufficiently large negative coupling to the measure term, or if several matter fields are 
coupled to the model. If this phase structure does indeed represent some behavior 
of the continuum it is plausible, based on the expectation of universality, that it 
should not depend on the details of the discretization and should also be present 
for degenerate triangulations. Thus corresponding simulations of modified models 
of three-dimensional simplicial gravity using degenerate triangulations could serve 
as a confirmation on the observed phase structure. This work is in progress. 

Although the result of my numerical investigation of maximally degenerate 
three-dimensional triangulations indicated that this ensemble does not lead to a 
well defined statistical model, I do not rule out the possibility that a less restricted 
ensemble of triangulations, than considered in this paper, could be used. The chal- 
lenge is to identify an appropriate easing of the restriction which defines a consistent, 
exponentially bounded, canonical ensemble. With the potential benefits provided by 
an even larger ensemble than restricted degenerate triangulations this is definitely 
worth investigating. 
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